RM3-TH/00-17 
ROMA-1297/00 



First Lattice Calculation of the 
Electromagnetic Operator Amplitude 

GLADIATOR 
The SPQcdR Collaboration 

Southampton-Paris- Rome 

D. Becirevic a , V. Lubicz 5 , G. Martinelli a and F. Mescia° 

a Dip. di Fisica, Univ. di Roma "La Sapienza" and INFN, Sezione di Roma, 
Piazzale Aldo Moro 2, 1-00185 Rome, Italy. 

b Dip. di Fisica, Univ. di Roma Tre and INFN, Sezione di Roma Tre, 
Via della Vasca Navale 84, 1-00146 Rome, Italy. 



Abstract: 

We present the first lattice calculation of the matrix element of the electromagnetic operator 
(tt°\Q+\K ), where Q+ = (Q d e/lQir 2 ) (s L a^F^d R + s R a^F^d L ). This matrix element 
plays an important role, since it contributes to enhance the CP violating part of the 
Kl 7c°e + e~ amplitude in supersymmetric extensions of the Standard Model. 



1 Introduction 



The origin of CP violation is one of the fundamental questions of particle physics and cos- 
mology which remains an open problem to date. The recent measurements of e'/e 0] have 
definitively established direct CP violation and ruled out superweak scenarios. Unfortu- 
nately, we are still far from a full quantitative description of the dynamics which generate 
the amount of CP violation observed in hadronic processes M. Given the large theoret- 
ical uncertainties affecting the calculation of e'/e, it is very useful to collect additional 
experimental information about CP violation in different processes. The most interesting 
ones are those for which CP violating effects are suppressed in the Standard Model (SM) 
and enhanced in its extensions. Among the processes which have been considered in the 
literature, we like to mention charge asymmetries in non-leptonic decays |3|, |J and CP 
asymmetries of hyperon decays f| . 

Good candidates to provide new large CP violating effects are the super symmetric 
extensions of the SM with generic flavour couplings and minimal particle content. In 
this framework, among the possible contributions, it has been recently recognized the 
importance of the electromagnetic and chromomagnetic operators (EMO and CMO) 

Qj = ^(s L a^F^d R ±rs R a^F lMU d L ) (1) 
Qf = T ^(s L a^t a G%d R ±s R a^t a Gl l/ d L ) (2) 

The same mechanism, the misalignment between quark and squark mass matrices, may 
indeed substantially increase their CP-odd contribution to physical processes. In previous 
studies, particular attention has been devoted to the CMO which, without conflict with the 
experimental determination of the K°-K° mixing amplitude, can account for the largest 
part of the measured e'/e |§-@. 

In this paper, we consider the CP violating contribution of the EMO to Kl — > 7r°e + e~. 
The master formula which has been used in the numerical calculation of the rate is 

B{Kl _ A+e - )EM0 = 5.3 x 10 -< ( |g^g^ ) 2 ^ {lmS+? . (3) 

The coupling 5 + is related to the splitting in the down-type squark mass matrix. The 
definitions of 6 + and y 7 can be found in sec. |2] and x gq = m^/m^ is the ratio of gluino and 
(average) squark mass squared. The numerical coefficient in eq. ([3]) is the appropriate one 
for the operator renormalized in MS at the scale \i = 2 GeV. 

Our main result is the first lattice calculation of the B parameter Bt (also defined in 
sec. 0), for which we obtain 

B^Qi = 2GeV) = 1.21 ± 0.09 ± 0.04+°;°£ , (4) 

where the first error is the statistical one, the second is the systematic error due to the 
uncertainty on the ratio of the EMO to the vector current matrix elements and the third 
is the error coming from the uncertainty on the renormalization of the magnetic operator. 
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From the experimental upper bound JTD 

B(K L 

by taking from eq. 



vre'e ) < 5.1 x KT 10 (5) 

1 and rrig = 500 GeV, and using eq. (f|), we obtain 

\Im5+\ < 1.0 x 10~ 3 (95% C.L.) (6) 

The remainder of the paper is organized as follows: in sec. H all the formulae necessary 
to derive eq. are presented; in sec. ^ we describe the lattice simulation and discuss the 
calculation of Bt\ sec. f| contains our conclusion. 



' x gq 



2 EMO Contribution to B(K L -> 7r°e + e ) 

In this section, we recall the main ingredients necessary to compute ~B(Kl — > 7r°e + e _ ) in 
SUSY and introduce all the quantities appearing in eq. ([|). We discuss separately the 
effective Hamiltonian and the calculation of the branching ratio. 



2.1 The effective Hamiltonian for the magnetic operators 

The supersymmetric contribution to the effective Hamiltonian, in the case of the magnetic 
operators, can be written as 

Hmo = C 7 + Gu)Q» + C 7 -(/i)Q» + C^Q^/i) + C;Qj)Q-Qm) , (7) 

where the operators are renormalized at the scale \i. The Wilson coefficients generated by 
gluino exchanges at the SUSY breaking scale are given by 

c>5) = ^ [(iiUi±(ifs) 



n 



na s (mg) 



12. 



Fo{x gq ) , 



(8) 



Here {5% R 



(M 



I mi denote the off-diagonal entries of the (down-type) squark 



mass matrix in the super-CKM basis fI2" |. The explicit expressions of -Fo(^) and Go(x) are: 

Ax (1 + Ax - 5x 2 + Ax In (x) + 2a; 2 In (x)) 



F (x) 
G (x) 



3(1 -x) 4 

x (22 - 20x - 2x 2 + 16x In (x) - x 2 In (x) + 9 In (x)) 
3(1 -x) 4 



(9) 
(10) 



with i*o(l) = 2/9 and Gq(1) = —5/18. In the following we will use the combinations 
<5± = ($lr)2i i (^lr)i2 = (^lr)2i ± (3rl)2i- These quantities are the natural couplings 
appearing at first order in any parity conserving (+) or parity violating (— ) observable. 
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In the (Q^,Q^) basis, using the leading order (LO) anomalous dimension matrix 



it is straightforward to derive 



7 



8 
3 

32 
3 



(11) 



Cf{n) = vC}(m~ 9 ), (12) 

where 

2/21 / , n \ 2/23 / / x \ 2/25 



' a s (m s ) \ ( a s (m t ) \ ( a s (m h ) \ 
s a s (m t ) J \a s (m b ) J \ a s (/j) J 



2.2 Calculation of B(K T , -> 7r°e+e 



(13) 



In order to compute the rate, besides the Wilson coefficient, we also need the matrix element 
of the operator Q+, which is usually expressed in term of a suitable B parameter || 

(^IQ+I^ ) = *#^«^t7W) • (14) 

With respect to the standard definition of B T , we have introduced the g 2 -dependent form 
factor 1Z T (q 2 ) (7^-^(0) = 1) to account for the dependence of the matrix element on the 
momentum transfer q = — p n - Note that, since we are using renormalized operators, 
Bt depends on both the renormalization scheme and scale. 

Neglecting lepton masses, and isospin breaking effects, we may write the following useful 
identity 

(7r°e+e-|Q+|K°) Q d gB T , + 

TMf) = A*m K f+tf) {ir ^\{^e){sru)\K^) , (15) 

where a is the electromagnetic coupling and f + (q 2 ) = f + (0)TZ + (q 2 ) is the vector current 
form factor defined as 



(n°\(sYu)\K+) = -L 



(16) 



Eq. ( |l5f) allows us to write B(Kl — > 7r°e + e )emo in terms of the K + semileptonic branching 
ratio 

B(K L AV) IW = 2 (£)' B(K + . (17) 



where we have followed the notation of ref. || by defining 

Qd 



ImA+y, = ™ ImC+ , (18) 



with 

A 9 ( X 9<l) = 5 + G (x gq ) 

~ F o{x 9 q) , oil ..-J 



y^m~ g ,x gq ) = — — ? =- rj 

Trig \J lLx F m K 



G (x 



+ 8 1 -rf 



(19) 



The definition of y 7 given above differs from that of ref. || by a factor —Bt, since it is 
preferable to separate the Wilson coefficient from the B parameter. 
We have introduced the "effective" B parameter Bt defined as 

Br=4h*T, r.f^MMC, (20) 

/ + (0) j^K-m.)* (i? 2 A (3/2)( ? 2)| 1?+ ( ? 2)|2 

where A(g 2 ) = {m 2 K + m 2 — q 2 ) 2 — Am^m 2 . T is the correction due to the different q 2 
dependence of the tensor and vector form factors. In the calculation of T, we have used 
the experimental determination of the q 2 dependence of the semileptonic decay rate Jl3| . 

TZ+(q 2 ) = 1 + ^q 2 , A+ = 0.0286 ±0.0022 (21) 



and the slope of the tensor form factor extracted from our lattice data (see sec. ^ 

TZ T (q 2 ) = 1 + ^q 2 , A T = 0.022 ± 0.001 . (22) 

Since the correcting factor is very close to one, T = 0.99, its effect is practically negligible 
for the value of the effective B parameter, Bt- For the same reason, the difference between 
our value of A + in eq. (^) and the experimental value of eq. ( pT|) is ininfluential to the 
determination of Bt- 

In order to compute B T , we also need B T and / + (0). To wit, using the data discussed 
in sec. [3], we have followed two different procedures: 



we have taken the value / + (0) = 0.978 from ref. [14|. This number was obtained 



by neglecting isospin breaking effects. Using our result for the B parameter, Bt = 
1.23 ± 0.09, and T = 0.99, we get B T = 1.25 ± 0.09; 

• we have computed on our data the ratio Bt/ f + (0) extrapolated to the physical meson 
masses, obtaining Bt/ f + {0) = 1-18 ± 0.09. With the same value of T as before, in 
this case we get Bt = 1.17 ± 0.09. 

The difference between the two different procedures is taken into account in the systematic 
error, so that we get 

B T = 1.21 ±0.09 ±0.04 . (23) 

To this result, we add a very generous estimate (±6%) of the systematic error (to be 
discussed in the next section) due to the renormalization of the EMO. In this way we 
arrive to the result quoted in eq. @) of the introduction. This result is consistent with 



previous estimates from refs. [15, 16 



We have now all the necessary elements for the calculation of B(Kl — > TT°e + e~)EMO- 
Using eq. (|TT|), the definitions ( jT9|) and the values given in table [I], we arrive to the master 
formula in eq. @ which has been used, together with Bt, to constrain Im<5 + . 
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Parameter 


Value and error 


G F 


1.16639 x l(T 5 GeV^ 


V us 


0.2196 ±0.0023 


m K 


0.498 GeV 


m n 


0.135 GeV 


r(K L ) 


(5.15 ±0.04) x 10" 8 s 


t(K + ) 


(1.2385 ± 0.0025) x 10 8 s 


B(K+ -> 7r°e+z/ e ) 


0.0485 ±0.0009 


oc s {M z ) 


0.119 ±0.003 



Table 1: Average and errors of the main parameters. When the error is negligible it has 
been omitted. 



3 Lattice calculation of the EMO matrix elements 

In this section, we describe the procedure followed to obtain, in our lattice simulation, 
the B parameter, Bt (/ + (0)), and the form factor, IZri^q 2 ) (JZ+(q 2 )), necessary to the 
computation of B(Kl — > ir°e + e~)EMO- Since the calculation of the form factors on the 
lattice has been discussed in several papers, see for example |HJ and references therein, we 



only give here the details which characterize the present study. 

All our lattice results have been obtained using a non-perturbatively improved ac- 
tion |Tj|. The relevant operators, namely the vector and tensor currents 



V,j, = Z v (l + b v ma) [qj^Q + ic v d v qa ^Q] , 

T^u = Z T (/j)(l + b T ma)[iqa^Q + Or{d v qj^Q - d^qjvQ)) , (24) 



are improved. In our study, the coefficient bx, computed at 0(a s ) |]19| , is evaluated by using 



boosted perturbation theory |2T| . Zt(h) was obtained in ref. |21] with the non-perturbative 
method of ref. ||22|| , in the RI-MOM scheme at the renormalization scale fi. Note that for 



the tensor current, at the NLO, the RI-MOM scheme in the Landau gauge coincides with 
the MS scheme. The other constants are taken from the most recent non-perturbative 



determinations |23|, 53 1 . In summary, we have used the following values 



Z v = 0.79 , b Y = 1.4 , c y = -0.09 , 

Z T (fjt = 2 GeV) = 0.87 ± 0.01 , b T = 1.2 , c T = 0.05 . (25) 

Since the perturbative value of Zj- is Zj-(fi = 2 GeV) = 0.934(5), we estimate that the 
systematic error due to the normalization of the lattice operator is less than 6 %. Also this 
effect is included in the evaluation of the systematic error. 

Our analysis is based on a sample of 91 independent quenched gauge- field configura- 
tions, generated at the lattice coupling constant (3 = 6/#q = 6.2, on the volume 24 3 x 48. 
We use the combination of values of the hopping parameter, K^j, which are given in tab. 0. 
In order to calibrate the lattice spacing, a, and to extrapolate the form factors to the phys- 
ical meson masses, we use the lattice plane method p5[ . We find a" 1 = 2.7 ± 0.1 GeV, in 
agreement with previous simulations. 
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3.1 Extraction of the form factor 



Using standard lattice techniques, we have extracted from suitable correlation functions 
the matrix element of the operator 

(ttV^I^ ) = i (f KV l - pM) ^ fT } q2) , (26) 

where we have introduced the form factor /t(<? 2 ) to make contact with the definition used 
in ref . [17 ] . From eq. (|14]) , one finds 



2f T (q 2 ) B T K T {q 2 ) ^ 



rriK + rn-K tuk 

so that /t(0) = B T for rriK = Tn % corresponding to degenerate quark masses, m s = m u d . 

Besides /t(<? 2 ), we have also considered the vector form factors f + (q 2 ) and f°(q 2 ) 
(/ + (0) = /°(0)) appearing in eq. fllfi|). /t(<? 2 ), f + (l 2 ) and f°(q 2 ) have been computed for 
Pk = at several values of the pion momentum, p w = {2n / L)(n X) n y , n z ) in lattice units. 
We have results for (n x ,n y ,n z ) = (0,0,0), (1,0,0), (2,0,0), (1,1,0), (1,1,1) and (1,2,0). 
The results for /r(<? 2 ) as a function of the dimensionless variable q 2 a 2 are shown in fig. [I]. 

At fixed quark masses, we fit the lattice form factors to the expressions 

fr(q 2 ) = M0) (l + «t q 2 a 2 ) , f + {q 2 ) = /+(0) (l + a + q 2 a 2 ) . (28) 



The slopes in eqs. (plD and (p2|) are given by A +i t = a +i Tm 2 ct 2 - For friq 2 ) we have only 
used the points corresponding to (1, 0, 0), (1, 1, 0), labeled as squares in fig. [TJ, because the 
quality of the signal in the other cases is rather poor, and gets worse as the quark masses 
decrease. For completeness, in the figure, we have also shown the other points which have 
not been considered for the fit. The results for fr{q 2 ) at different values of the "strange" 
and light quark masses, corresponding to the hopping parameters Kl and K\ respectively, 
are given in tab. 0. 



7r°|scr 



fJbV 



d\K°) 





K x 


Mo) 


0.1344 


0.1344 


0.95(4) 


0.1344 


0.1349 


0.86(4) 


0.1344 


0.1352 


0.90(5) 


0.1349 


0.1349 


0.84(5) 


0.1349 


0.1352 


0.81(5) 


0.1352 


0.1352 


0.80(6) 



Mp(K L , Ki 



Mp(K h Ki 



4.42(23) 
4.85(31) 
6.48(24) 
6.05(35) 
6.86(29) 
7.21(48) 



0.090(1) 
0.073(1) 
0.064(1) 
0.058(1) 
0.049(1) 
0.039(1) 



0.090(1) 
0.058(1) 
0.039(1) 
0.058(1) 
0.039(1) 
0.039(1) 



Table 2: /r(0) and the ctT parameter for different combinations of the mesons masses, Mp, 
in lattice units. 
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Figure 1: friQ 2 ) as a function of q 2 for Kl — Ki — 0.1344. Beside every lattice point we 
mark the components (n x , n y , n z ) of the three momentum given to the pion in lattice units 
(2n/La). 



3.2 Extrapolation to the physical point 

The values of /t(0) and in tab. ^| have been extrapolated to the physical point, corre- 
sponding to M P (K L , Ki) = rriKCt and M P (K h K{) = m n a, with the lattice-plane method of 
ref . [^5] . Two different formulae have been used: 

• we have ignored the SU(3) symmetry breaking corrections, due to the m s -m Uj d mass 
difference, by making a fit of the form 

y = C+LMj,(K L ,K l ) , (29) 

where y = /t(0) or ax- The results are 

/ T (0) = 0.77(6) C M0 ) = 0.67(8) L Mo) = 2.9(6) 
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(30) 

a T = 7.8(4) C aT = 9.9(6) L aT = -60(7) 
From figs. |] and we see that SU(3) breaking effects are small since all the points 



1.1 




Figure 2: /r(0) as a function of the squared pseudoscalar meson mass in lattice units. The 



full (dashed) line represents a fit of the lattice points to eq. ( ]2^ ) for all (degenerate) meson 
masses. 



are rather close to the straight lines of the fit. The two lines correspond either to a 
fit to all the points or to those corresponding to degenerate mesons (M P (K L , K{) = 
M P (K h Ki)) only. 

we have chosen for /r(0) and «y a fitting formula which accounts for SU(3) breaking 
effects 

y = C + LM 2 P {K L ,K l ) + LM 2 P {K h K l ) (31) 
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Figure 3: as a function of the squared pseudoscalar meson mass. The curves represent 
a fit of the lattice points to the eq. (pS)). 



In this case the results read 

/ T (0) = 0.78(6) C /r(0) = 0.67(8) L /t(0) = 3.1(7) L M0) = -0.14(72) 

(32) 

a T = 8.1(4) C aT = 9.5(7) L aT = -38(10) L aT = -20(4) 

We conclude that SU(3) breaking effects are very small for /t(0) and quite small for 
a T . In the following, for physical applications we will use the results in eq. (|32"|). 



Using the physical K° and it masses, from eq. (0) we obtain 

B T = 1.23 ±0.09, A T = 0.022 ±0.001 , (33) 
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at a renormalization scale \i = 2 GeV. If one needs the operator at another value of //, the 
slope will remain the same, whereas Bt scales according to the formula [26 



B T (fx 2 




4/(33-2 n/) r 



2 / 12411 - 126^ + 52/2^ a ^ 2 ) - q,(^) 
+ 9 V (33-2n/) 2 J 4n 



4/33 



1 + 



1379a s (/i 2 ) - a s (//i) 



2178 



7T 



Br(^i) 
(34) 



where we set n/ = 0, since we are working in the quenched approximation. In ref. |2T 
it was shown that even the one-loop evolution gives a very satisfactory description of the 
scale dependence of the matrix elements of the non-perturbatively renormalized EMO. 
We also present our result for the ratio B T j / + (0) extrapolated to the physical point 



Bn 



1.18 ±0.09 



(35) 



This number is lower than the result obtained by combining Bt = 1.23 ±0.09 and / + (0) 
0.978, namely Bt/ f + (0) = 1.26±0.09. This happens because, on our data, we get / + (0) 
1.04 ± 0.06 , which is 6% larger than the result of ref. 



14 . This is why it is important to 



have a direct determination of Bt/ f + (0)- systematic effects are expected to be smaller in 
the ratio. Moreover, the comparison between the two ways to compute Bt/ / + (0) allows 
us to evaluate the systematic uncertainty. Indeed the difference between the two way of 
determining Bt/ f + {0) is compatible in size with the uncertainty which can be estimated by 
extrapolating the data using different procedures. In this case, the different extrapolations 
give results varying by about 5%. Thus we take 5% as a measure of a further systematic 
uncertainty on the determination of the ratio -&r// + (0). From eqs. (|33D , using / + (0) = 
0.978, and from eq. fl35|, we obtain B T = 1.25 ± 0.09 and B T = 1.17 ± 0.09, respectively. 
Considering the difference as a systematic error, we end up with our final result, which 
was given already in eq. 



4 Conclusion 

We have presented the first lattice calculation of the matrix element (vr ^^!^ ). The 
operator is renormalized non-perturbatively in the RI-MOM scheme which, at the NLO, 
is equal to the MS scheme. Including the statistical and systematic uncertainties (except 
quenching), our final result has an error smaller than 10%. This allows, from the upper 
limit on B(Kl 7r°e + e~), to put more stringent bounds on squark-mass differences in 
supersymmetric extensions of the Standard Model. 
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